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NEUTRONICS ANALYSIS OF AN OPEN-CYCLE HIGH-IMPULSE 
GAS-CORE REACTOR CONCEPT 
by Charles L. Whitmarsh, Jr. 

Lewis Research Center 

SUMMARY 

In order to develop a procedure for calculating the critical fuel mass as a function of 
the propellant pressure, an analytical study was performed on a 196 600-newton- 
(44 200-lbf) thrust, 4400-second-specific-impulse open-cycle gas-core rocket reactor. 
Neutron transport calculations were made on a range of reactor configurations so that 
the required fuel mass and propellant pressure of a given reactor configuration could be 
calculated for various combinations of cavity diameter, reflector -moderator thickness, 
and amount of structural material included in the beryllium oxide (BeO) reflector - 
moderator. The ranges of variables calculated were 3. 048 to 4. 267 meters (10 to 14 ft) 
for cavity diameter, 0.457 to 0. 762 meter (1. 5 to 2. 5 ft) for reflector -moderator thick- 
ness, and up to 6 -volume -per cent separated molybdenum as a structural material in the 
reflector -moderator. The reactor model was a spherical cavity region consisting of a 
ball of plasma uranium -235 fuel suspended in a hydrogen propellant. The cavity was 
surrounded by a spherical shell regions of BeO cavity liner, feed hydrogen, BeO 
reflector -moderator, plenums for helium coolant, and a titanium pressure vessel. 

High -temperature (4160 to 22 400 K) hydrogen was observed to significantly affect 
the neutronics of this reactor. Negative reactivity worths of the order of 15 percent 
Ak/k were attributed to upscattering of neutrons by high-energy hydrogen atoms. Con- 
sequently, median fission energies in the range of 0. 3 to 0. 7 electron volt resulted from 
hardening of the spectrum of neutrons returning from the reflector to the core region. 

In order to increase the mass of uranium fuel contained in this reactor without in- 
creasing the uranium loss rate, it is necessary to operate at a high cavity pressure. 

The net reactivity effect (positive fuel and negative propellant) decreases with increased 
fuel loading. A limiting pressure was calculated which represents the pressure corre- 
sponding to the fuel loading at which this net reactivity effect becomes negative. 

Reactivity coefficients were calculated for use in analyses of the dynamic behavior 
of open-cycle gas-core reactors. 


INTRODUCTION 


The gas -core reactor concept is applicable to manned interplanetary travel because 
of its potential for high thrust -to -weight ratio and high specific impulse (refs. 1 to 3). 

Of interest in this analysis is the open-cycle coaxial -flow concept, in which a critical 
mass of plasma nuclear fuel is suspended in a reactor vessel by fluid dynamic forces of a 
hydrogen propellant gas flowing around it. This core is surrounded by a reflector - 
moderator to reduce fuel requirements. A schematic in figure 1 illustrates the reactor 
concept as used in this report. 



Figure!. - Schematic of open-cycle gas-core reactor. 


Although lower specific impulse (1500 to 2500 sec) concepts have been examined both 
analytically and experimentally for several years in the coaxial -flow gas -core reactor 
program at the Lewis Research Center, the concept of high-impulse engines (of the order 
of 3000 sec or higher) and the design changes required for them are relatively new. One 
major variation is the use of a high -temperature reflector material and a reflector cool- 
ing mechanism (ref. 4) to remove the approximately 7 percent of reactor power that is 
deposited through neutron and gamma ray absorption. At specific impulses less than 
3000 seconds, regenerative cooling was sufficient to remove the neutron and gamma ray 
heating. To further explore the feasibility of high -impulse gas -core rocket engines, a 
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design study was initiated to determine the major problem areas and to better define 
technical feasibility. This report presents the neutronic data generated during that study 
(ref. 5). 

Previous analyses by Hyland (ref. 6) and Ragsdale (ref. 7) indicated the importance 
of the dependent relation between fuel mass and hydrogen propellant pressure in estab- 
lishing nuclear criticality and fuel containment in a gas -core reactor. A procedure was 
therefore developed for calculating a fuel mass which would consider both of these con- 
straints. Thus the calculated fuel mass would provide nuclear criticality and be sur- 
rounded by a hydrogen propellant at a pressure sufficient to contain that fuel mass. The 
procedure would be useful as an engineering design tool since it would allow calculation 
of the critical mass - reactor configuration - propellant pressure relation for a given set 
of operating conditions. This feature has been lacking in previous gas-core analyses. 

Reported herein are the data from neutron transport calculations which were per - 
formed on a range of reactor configurations in order to develop this procedure for var - 
ious combinations of cavity diameter, reflector thickness, and amount of structural ma- 
terial. Also included are ancillary data obtained in this study which pertain to neutronics 
effects of high -temperature hydrogen and to reactivity coefficients useful in kinetics 
analysis. 

Ground rules established to facilitate the reactor calculations are as follows: 

(1) A thrust of 196 600 newtons (44 200 lbf) and a specific impulse of 4400 seconds 
were selected based on a mission analysis of a manned Mars trip. 

(2) A ratio of fuel diameter to cavity diameter of 0. 67 was selected based on exper- 
imental cold -flow tests of stagnant gas containment within a flowing gas (refs. 8 and 9). 

235 

(3) Uranium enriched in uranium -235 (U ) was used for fuel because of its avail- 
ability, since usage might exceed 10 000 kilograms per trip (ref. 4). An enrichmant of 
98 percent was based on the maximum available in the Atomic Energy Commission sep- 
aration work cost schedule (ref. 10). 


CALCULATIONAL MODEL AND PROCEDURE 
Configuration 

The basic components of fuel, propellant, reflector, and pressure shell were ar- 
ranged in a spherical model in such a way as to minimize critical mass while still pro- 
viding the required functions of fuel containment, cavity liner heat protection, and struc- 
tural heat removal. The design bases for this reactor configuration and its components 
are described in reference 5. The analytical model for the neutronics analysis is de- 
scribed in table I. The reactor cavity consists of a core region of plasma fuel, sur- 

238 

rounded by propellant hydrogen containing a seed material U , and bounded by a porous 
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TABLE I. - DIMENSIONS AND MATERIAL COMPOSITIONS OF SPHERICAL 


MODEL USED IN CRITICALITY CALCULATION 


Region 

— 

Average 

temperature, 

K 

Region 

thickness, 

cm 

Material 3 " 

Composition, 

o 

atm/cm 

Fuel 

50 000 

b 102 
c 122. 5 
d 143. 1 

98 U 235 
2 U 238 

Variable fuel loading 

Cavity hydrogen 
(propellant) 

10 600 

b 50. 4 
c 60. 3 
d 70. 2 

95 wt % H 
5 wt % U 238 

N H = 0. 290X10 21 
u238 yj 

N u = 0. 614x10 1 ‘ 

Cavity liner 

1 600 

1.27 

84.7 BeO 
b 10 H 

N Be = N° = e 0. 0614X10 24 
N H = 0.377X10 21 
T t238 14 

N u =0.8x10 

Feed hydrogen 

1 600 

1.27 

95 wt % H 
5 wt % U 238 

N H = 0. 00377X10 24 
N u = 0. 8 xlO io 

Reflector edge 1 


2.54 

94 BeO 

N Be = N° = 0.0682X10 24 

Reflector* 

1 600 

e 43. 18 
h 58. 42 
*73. 66 

84.2 BeO 
10. 5 He 

N Be = N° = 0.0635X10 24 
N He = 0. 147X10 21 

Hot plenum 
(coolant outlet) 

1 390 

1.27 

4. 7 BeO 
95 He 

N Be = N° = e 0. 0034X10 24 
N He = 0.00199X10 24 

Plenum divider 


1. 27 

94 BeO 

N Be = N° = 0.0682X10 24 

Cold plenum 
(coolant inlet) 

1 280 

1.27 

He 

N He = 0.00228X10 24 

Outer plenum wall 


1.27 

94 BeO 

N Be = N° = 0.0682X10 24 

Insulation 

800 

15. 24 

10 Zr 0 2 

N Zr = 0.00273X10 24 
N° = 0. 00558X10 24 

Pressure shell 

300 

15.24 

Ti 

N Tl = 0. 0566X10 24 


a Unless otherwise noted, all materials are in volume percent. 
b For 3.048-m- ( 10 -ft -) cavity -diameter configuration. 
c For 3.658-m- ( 12 -ft-) cavity -diameter configuration. 
d For 4.267-m- ( 14 -ft-) cavity -diameter configuration. 

Composition based on pressure of 90. 5 MN/m (400 atm) used as the reference 
case. 

f 

For some calculations the reflector contained heat -exchanger tubes made of 
separated Mo. 

g For 0.457-m- (1.5 -ft-) reflector -thickness configuration. 
h For 0.61-m- (2-ft-) reflector -thickness configuration. 
l For 0. 762 -m- (2. 5 -ft-) reflector -thickness configuration. 




wall liner. Relatively cold hydrogen from a feed system fills a region adjacent to the 
cavity liner and flows into the high -temperature cavity. The reflector region contains 
heat -exchanger tubes to provide heat removal by helium coolant gas. The tube material 
was assumed to be the same as the reflector material, BeO, though metallic tubes were 
also considered. A coolant volume equivalent to 10. 5 percent of the reflector volume 
was calculated to be sufficient to remove the heat generated by neutron and gamma ab- 
sorption in the reflector material. Coolant volume fraction was assumed constant for all 
models and the calculated value of 10. 5 percent was based on design limits of about 1. 4- 
MN/m^ (200-psi) pressure drop from inlet to outlet plenums and a 550 K (1000° R) tem- 
perature rise between tubes. Plenums for the outlet and inlet coolant gas were located 
outside the reflector in order to keep reflector material as close to the core as possible. 
A foamed zirconium oxide insulation layer was provided to help protect the pressure 
shell from overheating. The titanium pressure vessel provided containment for the pres- 
surized core materials. The rocket nozzle was neglected in the calculations because of 
its small size. (It was estimated to have a throat diameter of 3. 8 cm, compared to a 
cavity diameter of about 3.7 m. ) 

Based on a previous survey, cavity diameters of 3. 048 to 4. 267 meters (10 to 14 ft) 
and reflector thicknesses of 0.457 to 0. 762 meter (1. 5 to 2. 5 ft) appeared to be the range 
of interest. At lower diameters and reflector thicknesses the required containment pres- 

O 

sure would be excessive (>101.4 MN/m , or 1000 atm). Region dimensions and compo- 
sitions for the calculational model are itemized in table I. Thicknesses of all regions 
other than the cavity and the reflector are nominal and were selected as reasonable from 
the standpoints of fabrications and coolant pressure drop. 

Temperatures in the hydrogen and fuel cavity regions were based on a heat transfer 

analysis of a 3. 048-meter- (10-ft-) cavity -diameter reactor with a 0.457-meter (1. 5-ft) 

o 

reflector, operating at a pressure of 40. 5 MN/m (400 atm), and with a thrust of 196 600 
newtons (44 200 lbf) and a specific impulse of 4400 seconds. This temperature profile 
was assumed to be applicable to all configurations studied herein. The average hydrogen 
temperature of 10 600 K represents the temperature corresponding to the average hydro- 
gen atom in the cavity region and was obtained by volume averaging the product of tem- 
perature and density (see appendix A). At engine conditions, molecular hydrogen is 
largely dissociated. Theoretical data relating hydrogen density to temperature and pres- 
sure obtained by using the method described in reference 11 have been used to determine 
atomic densities of hydrogen for nuclear calculations. 

The model described in table I was selected for the neutronics analysis even though 

some of the materials and region temperatures are still to be specified. For example, 
238 

U was used as the propellant seed material primarily because vapor opacity data are 
available which are required in heat transfer calculations. Uranium-238 should be rep- 
resentative of the class of materials that could be used for this application. Also, the 
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temperatures of the cavity liner and feed hydrogen region were somewhat arbitrarily se- 
lected. Heat transfer analysis indicates that the cavity liner temperature is insensitive 
to propellant heating and therefore would probably be determined by the incoming hydro- 
gen. This inlet temperature of the hydrogen could be adjusted over a reasonable range, 
depending on what, if any, heat load is imposed on it (e. g. , nozzle and/or pressure ves- 
sel cooling) prior to entering the cavity liner. 


Procedure 

The calculational procedure was to select a reactor configuration and pressure and 
then calculate the critical mass. A fluid dynamics correlation was then used to deter- 
mine whether the calculated critical mass could be contained by the assumed propellant 
pressure. If not, a new value for pressure was selected and the procedure was repeated. 

Containment of fuel in an operating gas -core reactor must satisfy the following fluid 
dynamics criterion 


P 


0. 0038 


M 1.385^0. 383 t 0. 383 
F r *sp 



( 1 ) 


where 

2 

P cavity pressure, MN/m 

Mp fuel loading, kg 

F thrust, N 

I_„ specific impulse sec 

sp 

D c cavity diameter, m 

Vp/V c volume fraction of fuel in cavity 

This criterion is based on fluid dynamics and heat transfer conditions in the cavity. 
Equation ( 1) is based on a derivation in reference 7. The fuel mass must also meet the 
criterion of criticality, which for a parametric analysis is convenient to calculate in the 
form 
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M c = M(ref. ) 


% 

m 


— (H press.) 
k 


+ 


— (H temp. ) + 
k 


— (structural material)! (2) 

k J % Ak 

k 


where 

critical mass of reactor with reference conditions, kg 
reactivity effect of hydrogen pressure, percent Ak/k 


% — (structural material) 
k 

AM 
k 

The reference critical mass M(ref. ) is defined for calculational convenience as the crit- 
ical mass of a gas -core reactor with hydrogen conditions of 40. 5 MN/m (400 atm) pres- 
sure and 10 600 K temperature and with no structural material in the reflector. Any 
reactor with these conditions is hereinafter referred to as a reference reactor, regard- 
less of size. 

For a given configuration, M was calculated based on estimated values for hydro- 
gen pressure, hydrogen temperature, and structural material. Equation (1) was then 
utilized to evaluate whether the calculated M and the estimated hydrogen pressure P 
were consistent values for fuel containment. This procedure was iterated until 

m f =m c . 

Criticality calculations were performed with the neutron transport code TDSN 
(ref. 12) using the S 4 Pj 19 -group options in spherical geometry. The notation scheme 
for code options can be described as follows: S 4 represents the angular quadrature 
scheme of order 4 for solving the S n transport equations, P ^ represents cross sections 
obtained from a Boltzmann transport equation solution using the first two terms of a 
Legendre polynomial expansion, and 19 -group signifies the number of discrete energy 
groups for which the multigroup transport equation will be solved. The use of S^P ^ 19- 
group options in neutron transport calculations was determined to be adequate for the 
criticality analysis of a heavy -water (DgO) -reflected, gas -core (uranium hexafluoride) 
reactor critical experiment at room temperature (ref. 13). These code options were as- 
sumed to be applicable to the high-temperature reactor studied herein. Downscattering 


reactivity effect of using hydrogen temperature gradient in 
propellant region, percent Ak/k 

reactivity effect of addition of structural material to 
reflector, percent Ak/k 

reciprocal of fuel reactivity coefficient, kg/percent Ak/k 


M(ref) 

% Q (H prec«?) 
k 

% — (H temp) 
k 
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was allowed between all groups , whereas upscattering was limited to the seven thermal 
energy (E < 2. 38 eV) groups. 

Neutron cross sections were generated with the GAM -II (ref. 14) and GATHER -II 

(ref. 15) codes. Twelve fast group cross sections (E > 2. 38 eV) were flux weighted over 
235 

a U slowing -down spectrum for use in the core region and over a BeO spectrum for 
use in all regions outside the core. Seven thermal group cross sections (E < 2. 38 eV) 
were flux weighted over a BeO spectrum at 1600 K for use in all reactor regions. Also, 
the nuclei in each region were assumed to be at their appropriate temperature in an oper- 
ating engine. The use of high -temperature nuclei in the calculation (as opposed to no 
temperature correction) significantly increased the hydrogen scattering cross sections. 

At equivalent temperatures the velocity of protons (atomic hydrogen) is equal to that of 
neutrons. Since the cross section, or probability of reaction, is a strong function of the 
relative velocity between a proton and a neutron, this effect is particularly important in 
the cavity of a gas-core reactor, where protons can exist at temperatures of the order of 
30 000 K. The GATHER n code was used to calculate free gas scattering kernels for fuel 
and hydrogen regions. The resulting cross sections are then flux weighted over the spec- 
trum of interest, in this case BeO at 1600 K. 


DESIGN PROCEDURE 

Reactivity effects of the major variables of hydrogen pressure, hydrogen temper- 
ature, and percentage of structural material were calculated for a range of conditions. 
These reactivity effects were transformed to required changes in fuel loading by the use 
of fuel reactivity worths. The reference reactor critical mass was adjusted with the cal- 
culated change in fuel loading to obtain the critical mass for a given configuration. 


Reference Reactor Calculations 

To provide a basis for critical mass determinations, criticality calculations were 

performed for cavity diameters of 3.048, 3. 658, and 4. 267 meters (10, 12, and 14 ft) 

and for reflector thicknesses of 0.457, 0. 610, and 0. 762 meter (1. 5, 2.0, and 2. 5 ft). 

o 

Hydrogen conditions in the cavity were assumed to be 40. 5 MN/m (400 atm) and 
10 600 K. The fuel was uniformly distributed within the core region, which was assumed 
to be the central 30 percent of the cavity volume. No mixing was assumed to occur be- 
tween fuel and hydrogen in the cavity. 

Critical fuel density is shown in figure 2 to decrease with increasing cavity diameter 
(or core diameter) and with increasing reflector thickness. However, total fuel loading 
increases with increasing cavity diameter (fig. 3). 
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Fuel reactivity worths were determined for a wide range of fuel loadings with cavity 
diameter as a parameter (fig. 4). At the higher fuel loadings the analytical model con- 
tained structural material in the reflector region, and at lower fuel loadings the refer- 
ence reactor was used. However, the correlation appeared to be consistent, thus indi- 
cating fuel worth as a function of fuel loading to be essentially independent of material 
variations external to the reactor cavity. The general shape of the curves has been ob- 
tained experimentally on a D20-moderated gas-core critical mockup (ref. 16). 


Hydrogen Conditions 

The hydrogen in the reactor cavity region, because of its location between the fuel 
and the reflector, significantly affects the multiplication factor. Neutrons which have 
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been thermalized in the reflector and then returned to the fuel region are the major 
source of reactivity. However, while passing through the hydrogen, these neutrons are 
subject to both absorption and scattering collisions. Absorption reduces fission prob- 
ability by capturing the neutron before it can reach a fissionable atom. Scattering can af- 
fect criticality by changing the path of the neutron and by changing the energy of the neu- 
tron, thereby either reducing or enhancing its fission probability depending on the energy 
dependence of the fission cross section. If the hydrogen is hot (at energies greater than 
thermal), the tendency is to upscatter the neutrons, that is, transfer energy from protons 
(atomic hydrogen) to neutrons at lower energy levels. The near -Maxwellian distribution 
of neutrons emanating from the reflector is thus hardened (increased in energy) before it 
reaches the fuel. The effect is to reduce criticality by upscattering the neutrons into an 
energy range where the probability of capture (parasitic absorption) relative to fission in 
the fuel is greater than at lower energies. 

The actual hydrogen cavity pressure is a function of the fuel mass (eq. (1)). There- 
fore, to generalize the reference reactor critical masses, the reactivity effect was cal- 
culated by substituting hydrogen at 20. 3, 81. 1, and 121. 6 MN/m^ (200, 800, 1000 atm) 
into the cavity of reactor. The resulting change in reactivity is plotted in figure 5. For 
these calculations the hydrogen density as a function of pressure was obtained from fig- 
ure 14 (discussed in appendix B). Increasing slope (rate of reactivity change per unit 
pressure) with increasing cavity diameter results from the greater width of the cavity 
hydrogen region. (For all configurations the ratio of fuel diameter to cavity diameter 
was 0. 67. ) This greater width allows more proton -neutron collisions for a given pres- 



Figure 5. - Reactivity effect of hydrogen pressure in reference reactor configuration. 
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sure. Neutron absorption in cavity hydrogen as a function of pressure is shown in 
table II to vary almost proportionately with pressure. The reactivity effect is com- 
pounded, however, by the fact that this absorption selectively removes low-energy neu- 
trons, thereby hardening the neutron spectrum before it enters the fuel region. As a re- 

o 

suit the median fission energy, which was 0. 31 electron volt at 20. 3 MN/m (200 atm), 
increased to 0. 72 electron volt at 121. 6 MN/m 2 (1200 atm). Greater parasitic absorp- 
tion in the fuel as the spectrum hardens is illustrated by the ratio of neutron captures to 
fissions increasing from 0. 226 to 0. 234 as hydrogen pressure increases from 20. 3 to 
121. 6 MN/m 2 (200 to 1200 atm). 


TABLE n. - EFFECT OF CAVITY PRESSURE ON CORE FLUX 
SPECTRUM FOR REACTOR 


[Cavity diameter, 3. 658 m (12 ft); reflector thickness, 0.609 m 
(2 ft); hydrogen temperature, 10 600 K. ] 


Cavity pressure 

Median fission 
energy , 
eV 

Ratio of neutron 
captures to fis- 
sions in fueled 
region 

Absorption in 
cavity hydro- 
gen region per 
source neutron 

MN/m 2 

atm 

20. 3 

200 

0. 31 

0.226 

0.0075 

40. 5 

400 

.39 

.229 

.0156 

81. 1 

800 

.60 

.231 

.0325 

121. 6 

1200 

.72 

.234 

.0490 


All previous calculations in this report have contained the simplifying assumption of 
a uniform hydrogen cavity temperature. However, in an operating engine a gradient will 
exist between the relatively hot hydrogen-fuel interface and the cool hydrogen - cavity 
liner boundary. The calculated temperature distribution in the cavity hydrogen, shown in 
figure 6, decreases from 26 300 K at the fuel interface to 1600 K at the cavity liner. Any 
variation of hydrogen temperature with cavity diameters and pressures was assumed to 
be sufficiently small that the data in figure 6 were applied to all configurations considered 
herein. For calculational purposes the hydrogen region was represented by five spher- 
ical shells ranging in temperature from 4160 to 22 400 K (appendix A). 

The reactivity effect of zoning the hydrogen was -0. 25 percent Ak/k for a 3. 048- 
meter- ( 10-ft -) diameter, 0.61-meter- (2.0-ft-) thick reflector configuration (case 1) and 
-0. 50 percent Ak/k for a 3. 658-meter- ( 12 -ft-) diameter, 0. 61-meter- (2 -ft-) thick re- 
flector configuration (case 2). Spectral and criticality data for these two cases are com- 
pared to the uniform 10 600 K hydrogen calculations in table III. The negative reactivity 
effect is seen to come from increased neutron -capture -to-fission ratios in the fuel. 
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3X10 4 



Figure 6. - Caiculationai mode! representation of hydrogen tempera- 
ture distribution in cavity of gas-core reactor. Cavity diameter, 

3. 048 meters (10 ft); pressure, 40. 5MN/m ? (400 atm); thrust, 

1% 600 newtons (44 200 Ibf); specific impulse, 4400 seconds. 


TABLE m. - EFFECT OF HYDROGEN TEMPERATURE DISTRIBUTION ON CORE PROPERTIES 


Core property 

Cavity diameter, 3.658 m (12 ft); 
reflector thickness, 0.61 m (2 ft) 

Cavity diameter, 3.048 m (10 ft); 
reflector thickness, 0. 61 m (2 ft) 


Hydrogen zone at 
10 600 K average 
temperature a 

Hydrogen separated 
into five zones 
(fig. 6) 

Hydrogen zone at 
10 600 K average 
temper ature a 

Hydrogen separated 
into five zones 
(fig. 6) 

Multiplication factor 

1. 0039 

0.9988 

0.9995 

0.9970 

Median fission energy, 
eV 

0.39 

0.37 

0.36 

0.34 

Ratio of neutron captures 
to fissions in fueled 
region 

0.229 

0.236 

0.227 

0.234 

Absorptions in cavity 
hydrogen region per 
source neutron 

0.0156 

0.0165 

0.0127 

0.0133 

Reactivity worth of 
zoning, % Ak/k 

-0. 50 

-0.50 

-0. 25 

-0. 25 


a Temperature corresponding to average hydrogen density in cavity. 
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Since the five -zone hydrogen case has regions at temperatures both above and below the 
uniform 10 600 K calculation, both hardening and softening of the flux spectrum occurs. 
The result is to broaden the low-energy distribution (fig. 7). The net effect on median 
fission energy (table HI) was a reduction of 0. 02 electron volt. However, the greater im- 
portance of the upscattered neutrons on criticality is indicated by an increase in fuel- 
capture-to-fission ratio. 



I l l l l l ll l l l l 

3.7 1.4 . 5 .18.067 61 8.3 2.38 . 414 . 08 . 025 

Energy, MeV Energy, eV 


Figure 7. - Normalized flux spectra at hydrogen-fuel interface. 4> 3 J ^ ^ <p(u)du = 1. 

The small reactivity decrease due to hydrogen zoning indicates the adequacy of using 
the temperature of an average hydrogen atom in the propellant cavity region. This aver - 
age is obtained by volume averaging the product of density times temperature of the hy- 
drogen. Also, total hydrogen atoms are conserved. 


Structural Materials 

Scoping calculations indicated that the gas -core reactor was sufficiently sensitive to 
neutron absorber materials that no standard structural material with acceptable high- 
temperature properties could be used in practical quantities. Thus, the proposed struc- 
tural material was separated molybdenum (Mo) which would be used in the alloy TZM. 
The Mo would be isotopically enriched to >98 percent of Mo^ and Mo^ in order to re- 
duce the low-energy neutron absorption cross section to as low a value as possible based 
on the difficulties (and cost) associated with the separation process. The actual isotopic 
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distribution used was not included because of classification of production data. 

Results of criticality calculations on various reactor configurations indicated that the 
relative critical mass as a function of Mo content in the reflector -moderator was nearly 
independent of cavity diameter and of reflector thickness (table IV). Relative critical 
mass is defined as the ratio of the critical mass for a reactor configuration with a given 
Mo content to the critical mass for that same configuration with no Mo in the reflector - 
moderator. Initial calculations using 1. 7 -volume -per cent niobium in the reflector indi- 
cated that relative critical mass varied from 2. 59 to 2. 73 as reflector thickness varied 
from 0.457 meter (1. 5 ft) to 0. 762 meter (2. 5 ft). Later calculations using 3. 035- 
volume -per cent separated Mo showed that relative critical mass varied from 1. 64 to 1. 62 
as cavity diameter varied from 3. 658 meters (12 ft) to 4. 267 meters (14 ft). These 


TABLE IV. - STRUCTURAL MATERIAL EFFECT ON CRITICAL MASS AS 


FUNCTION OF CAVITY DIAMETER AND REFLECTOR THICKNESS 


Structural 

material 

Content 
in reflector, 
vol % 

Cavity diameter 

Reflector thickness 

Relative 

critical 

mass 

m 

ft 

m 

ft 

Niobium 

1.7 

3.658 

12 

0.4512 

1.5 

2.59 


1.7 

3.658 

12 

.61 

2.0 

2.68 


1.7 

3.658 

12 

.762 

2.5 

2.73 

Separated 

3.035 

3.658 

12 

0.61 

2.0 

1.64 

molybdemim a 

3.035 

4.267 

14 

.61 

2.0 

1.62 


6.07 

3.658 

12 

.61 

2.0 

2.68 


a Greater than 98 percent Mo®® and Mo 1 ®®. 


values were judged sufficiently constant to warrant a single correlation for relative crit- 
ical mass against volume percent material which would include all cavity diameters and 
reflector thicknesses in this study. Such a correlation for separated Mo is presented in 
figure 8 for a range to 6. 07 percent. 

The upper limit of 6. 07 percent was selected because this allows enough material to 
construct heat -exchanger tubes that could withstand 40. 5-MN/m^ (400-atm) external 
pressure (ref. 5). At higher material percentages the increase in critical mass appears 
to exceed the range of interest for design purposes. 
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Figure 8. - Effect on critical mass resulting from separated 
molybdenum being added homogeneously to reflector of ref- 
erence reactor configuration. 


Example Calculation 

To obtain the fuel mass for a specific configuration both the criteria stated in equa- 
tion (1) for containment conditions and in equation (2) for criticality must be satisfied. 
The reference configuration critical mass can be adjusted by the use of reactivity coeffi- 
cients developed in the previous sections to arrive at a fuel mass and hydrogen pressure 
consistent with the amount of structural material contained in the reflector. In order to 
utilize the structural material correlation in figure 8, equation (2) was rearranged as 
follows 


M c = M(ref. ) xR 


% — (H press. ) + % — (H temp. ) I 


Ak 

T 


( 3 ) 


where R equals the relative critical mass factor caused by structural material. 

Both hydrogen pressure and fuel mass are unknown initially. Once the cavity diam- 
eter, reflector thickness, and amount of structural material have been set, an iterative 
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procedure is required until consistent values are obtained. Equation (3) is solved for the 
appropriate values of M(ref. ), R, and H(temp. ) by using an estimated hydrogen pres- 
sure. The resulting is then substituted into equation (1) to obtain the hydrogen pres- 
sure P required for containment. If necessary, this procedure is iterated until P 
agrees with the estimated hydrogen pressure. 

An example problem follows for a 4. 267 -meter- (14 -ft-) cavity -diameter, 0. 61- 
meter- (2.0 -ft-) reflector -thickness configuration containing 1.9 -percent separated Mo 
in the reflector region: 

M(ref. ) = 63.6 kg (fig. 3) 

R = 1. 38 (fig. 8) 


% — (H temp. ) = -0. 65 (derived from table m) 
k 


The pressure is estimated to be 50. 7 MN/m' 


— (H press. ) = -3. 1 (fig. 5) 
k 


M = 63.6X1.38 + (0.65 + 3.1) %^_^L= 87.8 + 3.75 %^Ex_£M_ 

k <^Ak k pi Ak 

° k ° k 


(3) 


By using figure 3, the 4. 267 -meter curve is integrated from 87. 8 kilograms to some 

value of M c at which the area under the curve equals about 3. 75 percent Ak/k. Thus, 

M - 87. 8 + 3. 75 = 87. 8 + 13. 6 = 101. 4 kg. 

c 3.91 


p = o 0038 (101. 4) lj 385 ( 196 600)°’ 383 (4400)°- 383 

(4. 267) 4 ' 54 (0. 3) 1 * 51 

= 0. 0854 (101. 4) 1 - 385 =51.2 MN/m 2 

O 

Since 51. 2 MN/m exceeds the estimated pressure, the pressure is reestimated to be 
55. 8 MN/m 2 . 
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% — (H press. ) = -4. 6 (fig. 5) 
k 


M c = 87. 8 + 5. 25 |2i| = 87. 8 + 19. 9 = 107. 7 kg 
P = 0. 0854 (107. 7) 1 ’ 385 = 55. 8 MN/ra 2 

Since the calculated pressure equals the estimated pressure, the hydrogen pressure will 
contain the fuel mass. 

As a check on this method, this configuration was calculated directly. The calcu- 
lated k of 1. 0087, compared to k = 1. 0, indicated the adequacy of treating the para- 
metric effects of pressure, temperature, and structural material independently. 


Calculated Critical Mass 

Predicted reactor conditions of critical fuel mass and hydrogen pressure are quite 
sensitive to the accuracy of the calculated multiplication factor k, due to the low net re- 
activity worth of fuel over much of the design range. The significant effect of a few per- 
cent in k on the design conditions of a 3. 658-meter- ( 12 -ft -) diameter reactor with a 
0. 61-meter- (2. 0 -ft -) thick reflector is shown in figure 9. As k varies from 1 to 1. 05, 
critical mass varies from 44.3 to 58. 2 kilograms and hydrogen pressure from 32. 9 to 
48. 1 MN/m (325 to 475 atm). Constant -pressure lines are included in figure 9 to indi- 
cate the lesser effect on critical mass if the corresponding hydrogen pressure increase is 
not considered in the calculations. The relative change in critical mass increases as de- 
sign pressure increases because of the lower net reactivity effect discussed previously 
and tabulated in table V. 

A possible bias in the criticality calculations of a few percent could come from a 
number of sources. Critical experiments were necessarily performed at lower temper- 
atures, thereby omitting the reactivity effect of a hot hydrogen region between the fuel and 
the reflector. This analysis has shown that the resulting effect of high-energy hydrogen 
atoms on k is substantial due to spectrum hardening. Various potential sources of error 
in this study and their estimated effect on k are itemized in table VI. It is impossible to 
estimate with any confidence the net effect of these phenomena, although the potential 
exists (in either direction) for a major change in the numerical results of this study. 

Also neglected in this analysis is excess reactivity which would be required by a re- 
actor control system. This could be of the order of 4 or 5 percent Ak/k. If desired, 
excess reactivity could be included in the design procedure described herein by requiring 
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I 



Multiplication factor, k 


Figure 9. - Effect of calculated multiplication factor on critical 
mass in reference reactor configuration. Cavity diameter, 

3. 658 meters <12 ft); reflector thickness, 0.61 meter (2 ft). 


TABLE V. - NET REACTIVITY WORTH OF INCREMENTAL FUEL 
ADDITION TO GAS CORE REACTOR 


[Cavity diameter, 4.267 m (14 ft); thrust, 196 600 N (44 200 Ibf): 
specific impulse, 4400 sec. ] 










TABLE VI. - POSSIBLE UNCERTAINTIES IN CRITICALITY CALCULATIONS 


Phenomenon 

Comment 

Estimated change in 
reactivity, 

Ak/k 

Beryllium cross sections 

Possibly inaccurate based on previous reactor 
analyses; e.g. , fast critical experiments 
(ref. 18) and the Tungsten Water Moderated 
Reactor program at Lewis Research Center. 

-6 to +2 percent. 

Transport code options 
P order, energy groups, 
mesh spacing 

Original calculation technique established for 
a cold, heavy -water (DgO) -reflected reactor 
without atomic hydrogen. Upscattering caused 
by high -temperature hydrogen indicates that 
more energy groups are needed in the 
eipthermal range. 

Possibly a few percent 
either way. 

Reactivity coefficient 
interdependence 

Design procedure is based on independent ef- 
fects of hydrogen temperature, hydrogen pres- 
sure, and structural materials. The degree to 
which this is wrong will affect k. 

Possibly a few percent 
either way. 

Fission products 

Xenon- 135 buildup has been neglected. 

About -0. 8 percent. 

Material impurities 

Reactor -grade beryllium oxide can contain up 
to a few hundred ppm of impurities. 

Up to a few percent 
negative. 

Fuel shape 

Although a spherical shape was assumed, ex- 
perimental data indicate that it would probably 
be closer to a teardrop. 

Undetermined. 

Fuel -propellant boundary 
layer mixing 

Although a well-defined boundary was assumed, 
diffusion and mixing will probably occur. 

Undetermined. 

Fuel -diameter -to -cavity - 
diameter ratio 

Criticality calculations are sensitive to this 
1 value, which has not been well defined 
experimentally. 

i . 1 

Significant either way. 


that the critical mass and pressure correspond to k excess instead of k = 1. 

The calculational model used herein to represent a gas -core reactor is based on 
today's best guess of reactor conditions. However, it should be recognized that the cal- 
culated data in this study are particularly sensitive to the fuel -diameter -to -cavity - 
diameter ratio, which is probably known with less accuracy than any other reactor con- 
dition. If the fuel -diameter -to -cavity -diameter ratio were larger than 0. 67, reactivity 
would increase both because more fuel would be in the cavity and because less hydrogen 
would be in the cavity. Increased reactivity could be translated into a lower critical fuel 
loading for a given reactor configuration. Of course, the converse would be true if the 
fuel -diameter -to -cavity -diameter ratio were less than 0.67. 
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The calculational model was based on complete separation of fuel and propellant pri- 
marily because of a lack of fluid mechanics information. Based on qualitative experimen- 
tal results a more realistic analytical model would include propellant within the fuel re- 
gion. Some propellant could permeate the entire fuel region, and it is probable that a 
transition region of fuel and propellant exists between the separate flow regions. The re- 
activity effect of hydrogen inside the fuel region is positive. However, for a given fuel 
loss rate the net effect of hydrogen dilution on pressure is not known. Consequently, the 
phenomenon is recognized in this study but no attempt has been made to evaluate the ef- 
fect on required fuel loading. 


CORE CHARACTERISTICS 

Ancillary results were derived from the study required for development of the afore- 
mentioned design procedure relating critical mass to propellant pressure. These data 
are presented in this section. 


Limiting Pressure 

Inspection of equations (1) and (3) indicates that any addition of fuel to a subcritical 
configuration also requires the addition of hydrogen through increased pressure. Thus, 
the net reactivity effect (positive fuel reactivity and negative hydrogen reactivity) would 
be a better indication of real fuel worth to a design. As the fuel loading of a reactor in- 
creases, the reactivity worth of the fuel decreases rapidly; however, hydrogen worth is 
relatively constant. Consequently, the net reactivity effect decreases at a faster rate 
(table V). These data indicate that there is some pressure above which fuel is worth so 
little that the net reactivity effect of adding fuel is negative. Above this fuel loading (or 
propellant pressure) a critical engine cannot be obtained without the use of some control 
system. For example, above this limiting pressure the addition of fuel makes the reac- 
tor more subcritical because the negative worth of the hydrogen required to contain the 
fuel mass is greater than the positive worth of the fuel. Figure 10 shows the limiting 
pressure for 3.048-, 3.658-, and 4. 267 -meter- (10-, 12-, and 14 -ft -) diameter reactors 
to be 75, 69, and 62 MN/m 2 (720, 680, and 610 atm), respectively. Of course, these re- 
sults apply only to the range studied herein. It is conceivable that the reactivity worth 
curves in figure 10 could become positive at a much higher pressure, where the hydrogen 
region might act as a reflector. 

The preceding analysis was based on a reactor with no reactivity margin for a con- 
trol system. A similar analysis with excess reactivity would lead to larger limiting 
pressures. 
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Figure 10. - Net specific reactivity worth of fuel addition to reference re- 
actor configuration. Thrust, 196 600 newtons 144 200 Ibf); specific im- 
pulse, 4400 seconds. 


Hydrogen Upscattering and Temperature Effects 

As indicated previously, the upscattering of neutrons by high -temperature hydrogen 1 

atoms can be expected to reduce criticality. The magnitude of this effect was investi- 
gated by performing a criticality calculation using a special set of hydrogen cross sec- 
tions. This set included the same fast group cross sections as used in reference calcu- 
lations, but the thermal group microscopic cross sections were generated at 1600 K in- 
stead of 10 600 K. The macroscopic cross-section set was calculated from a hydrogen « 

density corresponding to 10 600 K. This procedure maintained the same absorption prob- 
ability (density dependence) as the reference calculations but isolated the upscattering ef- >| 

feet (temperature dependence). The results in table VII indicate that the upscattering ef- 
fect of 10 600 K hydrogen is worth about -15 percent Ak/k compared to 1600 K hydrogen. 

The negative worth results directly from the relatively greater capture probability in the 
fuel at higher temperatures (energies) and is indicated by the increase from 0. 218 to 
0. 229 of the ratio of neutron captures to fissions in the fuel region. The significance of 
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TABLE VH. - EFFECT OF HYDROGEN UPSCATTERING ON CRITICALITY 


OF GAS -CORE REACTOR 

[Cavity diameter, 3. 658 m (12 ft); reflector thickness, 0.61 m (2.0 ft); 
pressure, 40. 5 MN/m^ (400 atm). ] 


Treatment of hydrogen 

Multiplication 

factor 

Ratio of neutron 
captures to fissions 
in fuel 

Thermal group microscopic cross sections 
at 1600 K (density for macroscopic cross 
sections based on 10 600 K) 

1. 1878 

0.218 

10 600 K microscopic and macroscopic 
cross sections (reference calculations) 

1.0039 | 

.229 

Zoned regions (4 160 to 22 400 K microscopic 
cross sections with corresponding densities) 

.9988 

.236 


the upscattering worth of hot hydrogen to reactor design becomes apparent when com- 
pared to the previous section on the sensitivity of calculations to variations in k. 

The effect of hydrogen upscattering on criticality led to consideration of whether low- 
temperature hydrogen would offer a sufficiently lower fuel mass (and therefore lower 
weight) to warrant a lower specific impulse design. Calculations were performed with 
cavity hydrogen at 1600, 10 600, and 22 400 K and a constant pressure of 40. 5 MN/m^ 
(400 atm). The resulting data, itemized in table VUI, show that the lower critical mass 

TABLE VUI. - EFFECT OF HYDROGEN TEMPERATURE ON CORE 
PROPERTIES FOR GAS-CORE REACTOR 

[Cavity diameter, 3.658 m (12 ft); reflector thickness, 0.61m (2 ft); 
cavity pressure, 40.5 MN/m^ (400 atm). ] 


Core property 

Hydrogen average temperature, K 


1600 

10 600 

22 400 

Critical mass, kg 

-65.6 

46.5 

47.4 

Neutron absorptions in cavity hydrogen 
per source neutron 

-0.204 

0.0156 

0.0057 

Mean fission energy, eV 

-0.21 

0. 39 

0.30 

Ratio of neutron captures to fissions 
in fuel 

-0.210 

0.229 

0.254 
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occurred at 10 600 K and that only a slight weight penalty was incurred by raising the 
propellant temperatures to 22 400 K. However, in an operating engine this high- 
temperature weight penalty would be magnified. As propellant temperature is increased, 
the specific impulse increases and, as a result, the propellant pressure required to 
contain a given fuel mass increases (see eq. (1)). Greater hydrogen pressure would in- 
crease critical mass. The opposite argument can be made to show that the weight penalty 
in the 1600 K case is less than that calculated at constant pressure. To determine fuel 
mass as a function of propellant temperature would have required establishing a fuel- 
propellant pressure relationship at each temperature of interest, and that was considered 
beyond the scope of this report. 

At 1600 K, neutron absorption in the cavity hydrogen region was dominant; and at 
22 400 K, the upscattering of neutrons by the cavity hydrogen produced increased neutron 
capture in the fuel. The higher density of hydrogen at 1600 K (a factor of 13 higher than 
10 600 K) caused neutron absorption in the cavity hydrogen to increase from 1. 7 to 22. 2 
percent of the source neutrons. The lower density of hydrogen at 22 400 K caused a re- 
duction in neutron absorption (increased reactivity), but this effect was negated by the in- 
creased fuel -capture -to -fission ratio of upscattered neutrons. These data are not in- 
tended to indicate that the 10 600 K hydrogen is the best case but only to show the trend 
that neutron absorption becomes important at low temperatures and neutron energy level 
(or spectrum) becomes important at high temperatures. 

The decrease in median fission energy at a propellant temperature of 22 400 K ap- 
pears anomalous until one considers that several competing reactions are taking place - 
absorption, downscattering, and upscattering. Spectra at the fuel -hydrogen interface are 
plotted in figure 11 for propellant temperatures of 1600, 10 600, and 22 400 K. The high- 
energy peaks at about 2 MeV are neutrons from fissions in the core, and the low-energy 
peaks at 0. 2 to 2 electron volts are neutrons which have been thermalized in BeO and hy- 
drogen and returned to the core. Spectrum hardening due to temperature should shift the 
low-energy peaks toward higher energies as temperature increases. However, at con- 
stant pressure the hydrogen density decreases with increasing temperature, thereby al- 
lowing more low-energy neutrons to reach the core. When collisions do occur in higher 
temperature hydrogen, neutrons are upscattered to higher energies. This could explain 
the shape of the 22 400 K curve in which more low-energy and more high-energy neutrons 
occur than at 10 600 K. The net effect could be to lower mean fission energy. At 1600 K 
sufficient downscattering apparently occurs in neutron -proton collisions to more than 
compensate for absorption of low-energy neutrons. 

The hump in the 22 400 K spectrum in figure 11 at a lethargy of about 14 is not a real 
phenomenon but instead represents a calculation problem. The GATHER -II code used to 
generate microscopic cross sections has an upper energy cutoff of 2. 38 electron volts. 
Therefore, all upscattering reactions which cause the scattered neutron energy to exceed 
2. 38 electron volts are lumped into the next higher energy group (2. 38 to 8. 32 eV) instead 
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10 3.7 1.4 . 5 .18 . 067 61 8.3 2.38 . 414 . 08 . 025 

Energy, MeV Energy, eV 

Figure 11. - Flux spectra at fuel -hydrogen interface in reference reactor configuration with constant- 


temperature hydrogen in cavity at 40. 5 MN/m^ (400 atm). $ * / <p(ukiu * 1. 

*'-0.4 

of being distributed to their proper energy. For calculations at 10 600 K or lower, so 
few neutrons were scattered above 2. 38 electron volts that this procedure did not affect 
the transport calculation. At 22 400 K, however, a significant population of neutrons oc- 
curs in the 2. 38- to 8. 3 2 -electron -volt group. Better upscattering cross sections would 

tend to spread the neutrons out to higher energies , thereby smoothing out the curve in 
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figure 11. The effect on reactor calculations using U fuel would probably be to in- 
crease the multiplication factor because the average neutron energy would increase into 
an energy range where the ratio of capture to fission cross section decreases. (An ab- 
sorption resonance in the group above the cutoff energy results in a large capture -to - 
fission ratio in that energy range. ) The effect on reactor calculations with other fuels 
would depend on the behavior of their fission and absorption cross sections in the 2- to 
20-electron -volt energy range. In this study, however, the magnitude of any effect on k 
should be quite small because so few neutrons are affected in the reference configuration 
calculations (fig. 7). 


Reactivity Coefficients 


A number of reactivity coefficients were calculated for the gas-core reactor to sup- 
port a reactor dynamics study. These data, itemized in table IX, have been reduced to 
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TABLE IX. - CALCULATED REACTIVITY COEFFICIENTS FOR 


KINETICS ANALYSIS OF GAS-CORE REACTOR 


Property 

Calculated range 

Formula 


Limits 

Mean 



Lower 

Upper 

value 


Fuel mass, kg a 

42 

47 

44.5 

Ak/k = 0.38 (AM/M) 

Fuel temperature, K 

50 000 

60 000 

55 000 

Ak/k =0.001 (AT/T) 

Fuel radius, cm b 

122.6 

137.1 

129.85 

Ak/k = 0.21 (AR/R) 

Propellant mass, kg c 

8.7 

26.6 

17.6 

Ak/k = -0. 19 (AM/M) 

Propellant temperature, 
K 

1600 
1 10 600 

10 600 
16 200 

6100 
13 400 

AkA = -0. 11 (AT/T) 
AkA = -0.24 (AT/T) 

Moderator temperature, 
K 

1600 



1700 

1650 

AkA = -0.018 (AT/T) 


a Value is for 3.658-m- ( 12 -ft -) cavity -diameter configuration. Coeffi- 
cient varies with fuel loading and cavity diameter. Other values can 
be calculated from data in fig. 4. 

^Represents D v /D variation of 0.67 to 0.75. 

c r c o 

Corresponds to a pressure of 81.1 MN/rn (800 atm). 

the standard form Ak/k = F x Ap/p, where F is the coefficient and Ap/p is the rel- 
ative change of some property. Care must be taken, however, in interpreting the data 
because the coefficients have been specifically tailored for use in gas-core reactor equa- 
tions. Therefore, the calculation of each coefficient will be described in detail. In all 
cases the coefficient was determined from the reactivity difference between two static 
model calculations. 

Fuel mass. - The fuel density was varied between two calculations of the same 
configuration. 

Fuel temperature. - Only the effect of a temperature change on the microscopic 
cross sections of the fuel is considered here. In all calculations the low-energy neutron 
spectrum was assumed to be determined by the moderator temperature only. Thus, a 
fuel temperature change was represented by a nucleus temperature change in the calcula- 
tion of thermal cross sections. The result was a redistribution of the energy transfer 
cross sections. Uranium -23 5 was considered to be sufficiently like a 1/v absorber that 
the nucleus temperature had no effect on absorption cross sections (ref. 17). 

Fuel radius . - At a constant total fuel loading in the cavity, the fuel radius was var- 
ied. Since propellant volume was necessarily reduced by increasing the fuel radius, the 
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hydrogen density was varied to preserve the total number of hydrogen atoms in each 
calculation. 

Propellant mass. - At constant fuel loading and density, the propellant mass was var- 
ied by changing the hydrogen density in the cavity region. 

Propellant temperature. - The same procedure as for the fuel temperature coeffi- 
cient was used because the low-energy neutron spectrum in the propellant region was also 
assumed to be determined by the moderator temperature. 

Moderator temperature. - Again, only the effect of temperature change on micro- 
scopic cross sections was considered. However, a moderator temperature change was 
assumed to affect the neutron spectrum in all regions of the reactor. The BeO nucleus 
temperature change was also considered. The effect was to change low-energy absorp- 
tion and energy transfer cross sections in all reactor regions. 

As indicated, a basic assumption of the temperature effects on reactivity coefficients 
is that the low-energy neutron spectrum in all reactor regions is determined solely by the 
moderator temperature. Obviously this is not exact because of the upscattering effect 
observed in the hot hydrogen; but, judging from the fact that the fuel, feed hydrogen, and 
propellant hydrogen regions are all of the order of 1 or less mean free paths wide, the 
deviation is probably not large. 


Design Perturbations 

An interest in a different form of heat exchanger (pebble bed) led to the calculation of 
reactivity worth of a shell of material located at the inner edge of the reflector in the 
reference reactor. A 0. 625 -centimeter- (1/4 -in. -) thick shell of separated Mo was 
worth -19 percent Ak/k. This can be translated into critical mass increase by the use 
of fuel coefficients from figure 4 for the reference reactor configuration. For example, 
the reference critical mass of a 3. 058-meter- (10-ft-) diameter, 0.762-meter- (2.5-ft-) 
thick reflected configuration would increase from 29. 2 to 62 kilograms. The correspond- 
ing pressure would be 120 MN/m (1185 atm). The critical mass for this configuration 
as determined from equations ( 1) and (2) would be well beyond the range studied in this 
report. Therefore, any design configuration that would require metallic shells located 
between the reflector -moderator and the reactor cavity was not considered further. 

Also, this should give some indication of the reactivity penalty of using any structural 
material other than BeO (or other low -cross -section material such as beryllium or car- 
bon) to construct the cavity liner. 

Possible use of tungsten as the propellant seed material because of its high melting 

point led to calculation of its reactivity effect. For a reference reactor the reactivity 

238 

change was -0. 09 percent Ak/k when tungsten was substituted for U in the propellant 
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and feed hydrogen regions. Thus it would appear that the nuclear design is insensitive to 
the seed material since tungsten represents a relatively strong neutron absorber. 


CONCLUSIONS 

A procedure was developed for determining the critical mass as a function of propel- 
lant pressure of a coaxial -flow gas -core reactor for rocket propulsion. Data were gen- 
erated for a 196 600-newton- (44 200-lbf-) thrust, 4400-second specific impulse engine 
with a fuel -diameter -to -cavity -diameter ratio of 0. 67. Calculations were made for cav- 
ity diameters of 3.048, 3.658, and 4. 267 meters (10, 12, and 14 ft), reflector thickness 
from 0. 457 to 0. 762 meter (1. 5 to 2. 5 ft), a cavity temperature distribution of 1000 to 
27 000 K, and up to 6 -volume -percent separated Mo in the reflector for structural pur- 
poses. The following additional conclusions were drawn from these calculations: 

1. Consideration in the criticality calculations of the propellant pressure required 
for fuel containment significantly affects the predicted fuel mass in a gas -core reactor. 

2. Critical fuel loadings are sufficiently sensitive to parasitic absorption of neutrons 
that only small amounts (a few percent of the reflector -moderator volume) of separated 
molybdenum can be used as a structural material. Because of their higher neutron ab- 
sorption cross sections, other high -temperature structural materials could not be used. 

3. Spectrum hardening due to neutron upscattering in high -temperature hydrogen has 
a significant negative reactivity effect. Reactivity decreased 15 percent Ak/k when the 
hydrogen propellant temperature was increased from 1600 to 10 600 K at constant density. 

4. For a fixed set of reactor conditions of thrust, specific impulse, fuel -diameter - 

to -cavity -diameter ratio, and cavity diameter a limiting pressure (or fuel loading) exists 

above which a reactor cannot be made critical by increased fuel loading under operating 

o 

conditions. For the reference conditions this pressure is 73, 69, and 62 MN/m for cav- 
ity diameters of 3.048, 3.658, and 4. 267 meters, respectively. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 1, 1972, 

112-28. 
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APPENDIX A 


CALCULATIONAL TREATMENT OF HYDROGEN TEMPERATURE DISTRIBUTION 

Analytical representation of the hydrogen temperature distribution in the reactor 
cavity was accomplished by replacing the continuous distribution with five zones (spher- 
ical shells). Since temperature was assumed to be the more important variable (as op- 
posed to density), zone boundaries were determined by selecting approximately equal 
temperature intervals along the linear portion of figure 12 and closer intervals near the 
cavity liner. Corresponding hydrogen densities were also obtained from figure 12. The 
average temperature in each zone was obtained by 


_ y* T(r)p(r)r 2 dr 


f p(iV 


dr 


for each zone and the average density by 


_ f i p(r)r 2 dr 

P= 

/ r 2 dr 

where the integration limits are the zone radii. These data are itemized in table X. 
For larger diameter configurations the zone boundaries were scaled by conserving the 
volume fractions of each region. 


29 


.008 


.007 


.006 


.005 


^ * 


004 


.003 


.002 


.001 — 



Figure 12. - Hydrogen temperature and density in cavity of gas-core 
reactor. Cavity diameter, 3. 048 meter <10 ft); pressure, 40. 5MN/m 2 
(400 atm); thrust, 196 600 newtons (44 200 Ibf); specific impulse, 

4400 seconds. 


TABLE X. - FIVE -ZONE REPRESENTATION 


OF CAVITY HYDROGEN DISTRIBUTION 


Zone 

Average 

temperature, 

K 

Density, 

gem 3 

Boundary, m 

Inner 

Outer 

I 

22 400 

0.000178 

1.02 

1. 19 

II 

16 200 

.000296 

1. 19 

1.35 

m 

10 200 

.000531 

1.35 

1.47 

IV 

6 330 

.00111 

1.47 

1.51 

V 

4 160 

.00217 

1.51 

1.524 
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APPENDIX B 


HYDROGEN DATA 

At the temperatures and pressures of interest for a gas core reactor, hydrogen is in 
various stages of dissociation; and therefore the gas law is invalid for calculating gas 
density. For this report, densities for reference configuration calculations at 40. 5- 
MN/m (400 -atm) pressure are plotted in figure 13. Density variation with pressure at 
10 600 K is shown in figure 14. These data were calculated from procedures described 
in reference 11. 



Figure 13. - Hydrogen density at 40. 5 MN/m^ as function of 
temperature. 
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